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Abstract 

We investigate shadowing effects in deep-inelastic scattering from nuclei at 
small values x < 0.1 of the Bjorken variable. Unifying aspects of generalized 
vector meson dominance and color transparency we first develop a model for 
deep-inelastic scattering from free nucleons at small x. In application to nuclear 
targets we find that the coherent interaction of quark- ant iquark fluctuations with 
nucleons in a nucleus leads to the observed shadowing at x < 0.1. We compare 
our results with most of the recent data for a large variety of nuclei and examine in 
particular the Q 2 dependence of the shadowing effect. While the coherent inter- 
action of low mass vector mesons causes a major part of the shadowing observed 
in the Q 2 range of current experiments, the coherent scattering of continuum 
quark- ant iquark pairs is also important and guarantees a very weak overall Q 2 
dependence of the effect. We also discuss shadowing in deuterium and its im- 
plications for the quark flavor structure of nucleons. Finally we comment on 
shadowing effects in high-energy photon-nucleus reactions with real photons. 
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1 Introduction 



In recent years numerous experiments have been dedicated to high precision 
measurements of deep-inelastic lepton scattering from nuclei. Experiments at 
CERN |]-|J and Fermilab |5|-[| focus especially on the region of small values of the 
Bjorken variable x = Q 2 /2Mu, where Q 2 = —q 2 is the squared four-momentum 
transfer, v the energy transfer and M the nucleon mass. The data, taken over a 
wide kinematic range from 10" 5 < x < 0.65 and 0.01 GeV 2 < Q 2 < 100 GeV 2 , 
show a systematic reduction of the nuclear structure function F^(x,Q 2 ) with 
respect to A times the free nucleon structure function F^(x, Q 2 ) at x < 0.1. 

This so-called shadowing effect has prompted a fair amount of theoretical 
actvity (for recent reviews see e.g. ||). Some of the existing work focuses on 
an infinite momentum frame description of the scattering process (see e.g. ||10|| ). 
The driving mechanism in these models is given by quark and gluon annihilation 
processes at high parton densities, which are described using perturbative tech- 
niques. Although these methods allow one to address the Q 2 dependence of the 
shadowing effect, its x dependence is not accessible to perturbation theory and 
therefore subject to parametrization. 

Another class of models considers the deep-inelastic scattering process in the 
laboratory frame where the target is at rest [|IT| -|20]j . In this frame the interaction 
at small values of x proceeds via hadronic components present in the wave func- 
tion of the exchanged virtual photon. The coherence length of these hadronic 
fluctuations is typically of order 1/Mx and exceeds, for x < 0.1, the average 
nucleon-nucleon distance in nuclei. Hence for small x the hadronic fluctuations 
will interact coherently with several nucleons inside the target nucleus. Shadow- 
ing is then caused by destructive interference of multiple scattering amplitudes 
which describe the passage of these fluctuations through the nucleus. 

In this paper we present a laboratory frame description of deep-inelastic scat- 
tering at small x, based on ideas which unify the generalized vector meson dom- 
inance picture [^] with the concept of color transparency (for recent reviews 
see p2|). At small momentum transfers, Q 2 < 1 GeV 2 , the hadronic components 



of the absorbed virtual photon are formed by strongly correlated quark-antiquark 
pairs, most prominently by the low mass vector mesons p, u and <fi. At larger 
Q 2 quark-antiquark pairs from the so-called gg-continuum become increasingly 
important. We include both strongly correlated and continuum gg-fluctuations in 
terms of the measured photon spectral function. While some empirical informa- 
tion is available about the interaction of low mass vector mesons with nucleons 
and nuclei, the interaction properties of continuum quark-antiquark pairs are 
scarcely known. To fill this gap we use color transparency as a guiding prin- 
ciple, i.e. the cross section of color singlet quark-antiquark pairs is assumed to 
be proportional to their transverse size. With these ingredients we obtain a 
good description of both nucleon and nuclear structure functions, F^(x, Q 2 ) and 
F^(x,Q 2 ), at small x. 
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Although some of the ideas mentioned above are common to several recent 
models of deep-inelastic scattering at small x, little effort has been directed to- 
wards a quantitative comparison with the now available large amount of experi- 
mental data. We confront our model with most of the recent data. In particular 
we examine the Q 2 dependence of the shadowing effect - an intensely discussed 
issue. We find that while a major part of the shadowing seen in current experi- 
ments is caused by the coherent interaction of low mass vector mesons, coherent 
scattering of continuum qq pairs is also important and guarantees a weak Q 2 
dependence of the shadowing effect. 

We will also briefly discuss shadowing effects in high-energy photon-nucleus 
reactions with real photons (i.e. in the limit Q 2 -> 0). 

The plan of this paper is as follows: In Section 2 we introduce the space-time 
picture of deep-inelastic scattering in the laboratory frame. First we develop a 
model for deep-inelastic scattering from free nucleons in Section 3. Its extension 
to nuclear targets is described in Section 4. We discuss shadowing in nuclei with 
intermediate and large masses as well as its implications for a deuterium target. In 
Section 5 we apply our model to high energy photon-nucleus scattering. Finally, 
Section 6 contains a summary and conclusions. 



2 Lab frame picture of deep-inelastic scattering 

It is common to discuss deep-inelastic lepton scattering on free nucleons in a 
frame where the target moves with a large momentum, |p| — > oo. In this infi- 
nite momentum frame the parton model allows the interpretation of measured 
structure functions as momentum distributions of quarks and antiquarks in the 
target. There is, however, no reliable approach for dealing with nuclear systems 
in this frame. Consequently the preferable frame of reference for an investigation 
of nuclear effects in deep-inelastic scattering is rather the laboratory system in 
which the target is at rest. Well established knowledge about the structure and 
geometry of nuclear targets can then be used. 

Consider therefore a description of deep-inelastic lepton scattering from nu- 
cleons or nuclei in the laboratory frame. Here the basic photon-nucleon interac- 
tion process involves the time orderings shown in Figs. |I](a) and|l|(b): the photon 
either hits a quark (or antiquark) in the target which picks up the large energy 
and momentum transfer, or the photon converts into a quark-antiquark pair that 
subsequently interacts with the target. 

For small x the pair production process (b) dominates the scattering ampli- 



tude, as can be shown in time-ordered perturbation theory []I3"| : The amplitudes 
A a and Ab of processes (a) and (b) are roughly proportional to the inverse of 
their corresponding energy denominators AE a and AEb. For large energy trans- 
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fers v 3> M these are: 



AE a = E^-E^tx) , (1) 



y 2 , I(pI) + Q 2 

AE b = E b (t 2 ) - E h {t x ) « , (2) 

/ -A 1 / 2 

where (pi) is the average quark momentum in a nucleon and /x is the invari- 



ant mass of the quark-antiquark pair. We then obtain for the ratio of these 
amplitudes: 




(3) 



As we will argue later, the main contribution to process (b) comes from quark- 
antiquark pairs with a squared mass /i 2 ~ Q 2 . The ratio in Eq. (§) is evidently 
small compared to unity for x <C 0.1. Hence pair production, Fig. |l](b), will be the 
dominant lab frame process in the small-x region. For the following discussion 
of deep-inelastic scattering at small x we will therefore consider the dominant 
process (b) only (although it should of course be noted that in principle only the 
sum of (a) and (b) is Lorentz and gauge invariant). 

What are the implications of this picture for deep-inelastic scattering from 
nuclear targets? The coherence length A of a photon-induced hadronic fluctuation 
with mass /i is given by the inverse of the energy denominator (|2]) : 

A ~ J- = 2U ^ (4) 
AE b n 2 + Q 2 2xM { } 

For x < 0.05 this coherence length exceeds the average distance between nucleons 
in nuclei, d ~ 1.8 fm. Then coherent multiple scattering on several nucleons in 
the target can occur, leading to nuclear shadowing. 

For larger x, the coherence length of the intermediate qq state is small, A < d. 
At the same time the process in Fig. [TJ(a) becomes prominent, i.e. the virtual 
photon is absorbed directly by a quark or antiquark in the target. In the range 
of intermediate and large x, say x > 0.2, the virtual photon therefore interacts 
incoherently with nucleons bound in the nuclear target. 

3 Deep-inelastic scattering from free nucleons at small x 

Before discussing effects in the scattering from nuclei we have to develop a model 
of deep-inelastic scattering from free nucleons. The free nucleon structure func- 
tion F^(x,Q 2 ), defined as the average of the proton and the neutron structure 
function, can be written in terms of the virtual photon- nucleon cross section <t 7 * n'- 

F»(x,Q 2 ) = ±^-^a rN . (5) 



1 + % 4tt 2 



o ( , 
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In the limit v 3> Q and x < 0.1 that we are concerned with, this simplifies to 

F 2 N (x,Q 2 ) = -^a rN (6) 

47T Q. em 

(we use a em = e 2 /47r = 1/137). We will consider only contributions from trans- 
versely polarized photons, as they constitute the dominant part of the cross sec- 
tion: <t 7 *at ~ tf*N ( see HI] f° r an experimental analysis of <Jl/<Jt)- 

As discussed above, the virtual photon interacts with the nucleon by first 
converting into a qq pair which then propagates, forming a hadronic intermediate 
state that interacts strongly with the nucleon. This is expressed in the following 
spectral ansatz for the structure function |TT|,|14],pO|, 24], P^1 valid at x < 0.1: 

F 2 N (x, Q 2 ) = ^ f <W). (7) 

7T J Ami (/U 2 + Q 2 ) 

Here H(p 2 ) is the spectrum of hadronic fluctuations with mass p which is related 
to the measured cross section for e + e~ — > hadrons by 

n/ \ 1 ^"e+e~— >hadrons (^) /n \ 

S 



12tt 2 cr e +e-^ At + At -(s) ' 

Note that the effective hadron-nucleon cross section (JhNi.li 2 ) in Eq. (0) is an 
average including all contributions with a given invariant mass fi. The factor (/x 2 + 
Q 2 )~ 2 in Eq. ([/]) comes from the propagators of the hadronic intermediate states. 
One should of course note that Eq. (|7|) cannot be expected to follow directly from 
perturbative QCD. In particular, the effective cross section ^^(/x 2 ) incorporates 
non-perturbative physics characteristic of the small-x region. However, as will be 
shown, Eq. (0) does have the proper logarithmic behavior at large Q 2 . 

The structure function F 2 (x, Q 2 ) in Eq. (0) is dominated by contributions 
from intermediate states with an invariant mass /j, 2 ~ Q 2 . As a consequence for 
small momentum transfer, Q 2 < 1 GeV 2 , the low mass vector mesons p, uj and 
cf) are of major importance. They represent strongly correlated quark- ant iquark 
pairs which contribute the term 




n VMi V)= E ( — 1 W-O (9) 

v=p,uj,4> 

to the photon spectral function. Here my are the vector meson masses and 
gy l the corresponding coupling constants (see Table [1]). Eq. @ represents 
the traditional Vector Meson Dominance (VMD) model. For large Q 2 the heavy 
vector mesons J/ip and ip' also contribute, and we take them into account as well. 
Altogether, vector mesons give a contribution to the nucleon structure function 
of the form 

Ff™(x,Q 2 ) = ^E(— V ( Z^TT^V °vn- (10) 
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of the low mass vector mesons dominate the nucleon structure function FJ? and 
lead to the scale breaking behavior F%(x, Q 2 ) ~ Q 2 for Q 2 — > 0. 



Here <jvn are the vector meson-nucleon cross sections. They can be determined 
in real and virtual photoproduction experiments (see Table |l|). It should be 
mentioned that their exact value may in principle depend on the kinematics of 
the experiment (see e.g. [f29fl). 

For small values of x and Q 2 (x < 0.1 and Q 2 < 1 GeV 2 ) the interactions 

2 (x, Q 2 )~Q 2 for Q 2 ^0. 
For larger values of the momentum transfer, i.e. Q 2 > m 2 w 1 GeV 2 , the nu- 
cleon structure function F^ is governed by the interaction of quark-antiquark 
pairs with mass /i 2 ~ Q 2 > 1 GeV 2 . Apart from the narrow charmonium 
and upsilon resonances, these quark pairs form the so called qq continuum. 
In the annihilation of e + e~ into hadrons they are responsible for the approxi- 
mately constant behavior of the cross section ratio at large timelike momenta, 
a e + e -^ hadrons /<7 e + e -^ A1 + A1 - ~ 3J2f e /> where we sum over the fractional charges e/ 
of all quark flavors which are energetically accessible. 

To calculate the contribution of the continuum quark-antiquark fluctuations 
to the nucleon structure function we need to know their effective interaction cross 
section. Since the qq fluctuations of the photon are color singlets, we assume 
their cross sections to scale with their transverse size p (i.e. their size in a plane 



perpendicular to their momentum) as a ~ p 2 [ 16fl . Investigating the geometry 
of the dissociation of a photon into an uncorrelated qq pair more closely (see 
e.g. pO| ) one obtains the following approximate expression for p 2 : 

2 1 V 

fa . (11) 

1 a(l-a) (/i 2 + Q 2 ) 2 1 ; 

Here a is the fraction of the light-cone momentum carried by the quark: if k q 
is the quark momentum and q = (go, 0j_, 53) the photon momentum, one has 
& — {k q o + £^3) /(go + 53)- (Correspondingly, 1 — a is the light-cone momentum 
fraction carried by the antiquark.) Of course Eq. flTT| ) is a reasonable estimate for 
the size of the qq fluctuation only as long as the distance p of the pair is smaller 
than a typical confinement scale of about 1 fm. If the distance p increases further, 
strong interactions between the quark and antiquark will limit the transverse size 
of the qq fluctuation, thus leading to a saturation of p. Having this in mind, we 
choose for the effective cross section of continuum quark-antiquark pairs: 

f Rl 

a hN (p 2 ,a) = K ■ p 2 = K -mini 1 4/i 2 (12) 

1 a(l-a) (/i 2 + Q 2 ) 2 ' 

with a constant K to be determined. Here we have introduced a maximum ra- 
dius R c which should be in the range of the confinement scale. As it is clear 
from our discussion above, the cross section <ThN(p 2 , «) depends not only on 
the invariant mass p of the qq pair, but also on a, i.e. on the way the pho- 
ton momentum is split between the quark and antiquark. From Eq. (O) we 
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observe that the average interaction cross section of qq pairs with mass p is 
cr/j7v(/i 2 ) = Jq 1 da cr(/i 2 , a) ~ 1/p 2 (ignoring terms ~ log/i 2 ), which is the behav- 
ior necessary for scaling Hl4| , p!6| , p0| , |25|| . 

If we now take into account both the vector mesons and the quark-antiquark 
continuum, we obtain from Eqs. ([7], [10], [L2|) the following expression for the nu- 
cleon structure function: 

F 2 N (x,Q 2 )- 




, Q 2 r a 2 f 1 A /i 2 n cont -(/i 2 ) 2 

vr -W Jo (p 2 + Q 2 ) 

valid at small Bjorken x. Here n cont ' = II — n VMD is the continuum part of the 
photon spectral function which starts at Pq ~ m 2 . While the vector meson part 
vanishes as l/Q 2 for large Q 2 , the qq continuum contribution to the structure 
function displays logarithmic scaling behavior: 



F 2 N {x, Q 2 ) ~ In (R 2 C Q 2 ) for Q 2 > 1 GeV 2 . (14) 

We now compare our result for the free nucleon structure function with 
recent data of the New Muon Collaboration 0. We include in Eq. (O) all vector 
mesons p, u, <p, J /if} and if)' . Their masses, coupling constants and cross sections 
are summarized in Table |l|. 

The effective cross section cjvn (i-e. the forward scattering amplitude) may 
depend on the momentum and energy transfer variables Q 2 and v. The relevant 
range in Q 2 is, however, restricted by the fact that the vector mesons contribute 
mainly in the region Q 2 w my. On the other hand experimental constraints |].|2| 
put bounds on the accessible values of v. We therefore chose the cross sections 
to be approximately constant, setting o p n = 22 mb and fixing the other cross 
sections to scale like ctvn ~ V m y- 

The constant K in Eq. fll2|) is fixed at K = 1.7 together with R c = 1.3 fm. 
This corresponds to a maximum value of about 29 mb for the effective cross 
section of a gg-pair interacting with a nucleon. 

From Fig. ^ one can see that our model reproduces the measured nucleon 
structure function at small x quite well. We want to emphasize again the impor- 
tance of vector mesons at small values of Q 2 . In detail we find that at Q 2 = 1 GeV 2 
almost half of F 2 at x — 0.01 is due to vector mesons. At Q 2 = 10 GeV 2 they 
still contribute around 15%. 

4 Deep-inelastic scattering from nuclei at small x 

Just like the scattering from free nucleons, scattering from nuclear targets at 
small values of x proceeds via the interaction of hadronic components present in 
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the spectral function of the exchanged photon. For x < 0.1 the nuclear structure 
function can therefore be written in a way analogous to F^ in Eq. (|13D: 



F 2 A (x,Q 2 




T a hA^ 2 ,a)- (15) 



0^ 2 + Q 



We have just replaced the hadron-nucleon cross sections ahN in Eq- (0) by the 
corresponding hadron-nucleus cross sections a^A- 

As mentioned in Section |2], for x < 0.05 the coherence length A = 2v/ {Q 2 +/J?) 
of the interacting hadronic fluctuation exceeds the average internucleon distance 
in nuclei. Consequently the intermediate hadronic system can scatter coherently 
from several nucleons in the target. Interference between the multiple scattering 
amplitudes causes a reduction of the hadron-nucleus cross sections compared to 
the naive result of just A times the respective hadron-nucleon cross section and 
thus leads to shadowing. 

These effects are described by the Glauber-Gribov multiple scattering formal- 
ism |yj which we will now summarize briefly. 



4.1 Glauber-Gribov multiple scattering theory 

Let us consider high energy forward scattering of a hadronic fluctuation h with 
four-momentum q = (u, 0j_, ^v' 1 + Q' 2 ) and mass fi on a nucleus. In the labora- 
tory frame the target momentum is P = (A(M — £ ), 0), where A is the nuclear 
mass number and £ the binding energy per nucleon. The scattering amplitude 
AhA for this process can be written as the sum Aha — En=i Ah over multiple 
scattering terms A h n \ each of which describes the projectile interacting consecu- 
tively with n nucleons in the target (see Fig. §): 



(n) 



A\ 



(A-n): 
x f d 2 b 



n-l 

n 

i=l 



dL 



(2vr)2(M - £) 



dzi / dz 2 

J Z\ 

n-l 



dZr, 



(z/, . . . q iz . . .) 



xp„(b, z\ . . .z n ) Y[ [c 

i=l 



ili(zi-z i+1 ) 



(16) 



Here V^ n \u, . . . qi Z . . .) describes the interaction of the hadronic projectile with 
n nucleons. For large projectile energies u, it is assumed that V depends only 
on v and qi Z , the longitudinal momenta transferred to the interacting nucleons. 
Since we consider forward scattering only, we have Yh=i Qiz — 0. Furthermore 
the integration variables U are defined as U = Y^ l j=i<ljz, such that |q| — k is the 
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longitudinal momentum of the projectile after its interaction with the ith nucleon. 
The nucleon distribution in Eq. ([!(]) is given by the n-particle density 



p n (h;z 1 ...z % 



A 

n 

j=n+l 



d^x, 



5 3 (X r 



Sff(b,z 1 ',...;b,z n ',x ri +i...x A ) 



(17) 



where \1>(. . . Xj . . .) is the coordinate-space wave function of the nucleus. Its center 
of mass X cm = J2t=i x « i s fixed at the origin. Since the high energy scattering 
process occurs at a fixed impact parameter b, the active nucleons enter p n at 
coordinates x; = (b, z^) for i = 1 . . . n. 

As a next step the amplitude is expanded in hadronic eigenstates. Let us 
denote the complete set of states that can be reached after the interaction with 
the ith nucleon by {h{\ and write the corresponding invariant masses as m^. If 
the conversion from state hi into state hi + \ in the interaction with the (i + l)th 
nucleon is described by the transition amplitude fh,h i+1 , the expression for 
becomes: 



iV, 



{n) 



hh\ " 



X 



|q|-i 2 ) s 



(|q| -in-if 



(|q| ~k) 2 -m 2 hi +ie 
— -' l fh 2 h 3 x • • • 

ifhn-ih- 



ifhihi 



m 



1>2 



ie 



mi , + ie 



We can now perform the integration over the variables U. We note that the 
exponential factors in Eq. ([16]) require that the integration contour be closed in 
the lower plane. Picking up the poles, the longitudinal momentum transfer gets 
fixed at 



h 



|q 



v — mi 



Q 2 + ml 



(19) 



For intermediate and heavy nuclei we may in good approximation consider 
only elastic rescattering of the incoming hadronic state h from the nucleons inside 
the target. Contributions of inelastically produced states to multiple scattering 
were investigated by Murthy et al. [^2| and Nikolaev |33] for high energy hadron- 
nucleus scattering processes. They found such contributions to be small, though 
rising logarithmically with the projectile energy v. For example at v ~ 100 GeV 
inelastic terms account typically for ~ 5% of the total hadron-nucleus cross sec- 
tions under consideration. 

In the so-called "diagonal approximation", i.e. neglecting inelastic intermedi- 
ate states, the amplitude V^ 1 reduces to: 



(n) 



h diag. 



n 



hh- 



q| 



UY - /i 2 + ie. 



i/i 



hh i 



(20) 
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In this case the longitudinal momentum transfer is fixed just at the inverse of the 
coherence length A (B) of the hadronic projectile: 



i = |q| 



2v 



(21) 



Summing over all multiple scattering terms A£ and neglecting the binding 
energy ^cMwe find for the hadron-nucleus forward scattering amplitude 



iA, 



hA 



E 



.4! 



{{(A-n)\ 



(4Mv)- n+1 (if hh ) 



X / d 2 b / d Zl / dz 



dZ r . 



xp n (b;z 1 . . .z n )exp i 



. Zi Zn 



A 



(22) 



With the assumption that the hadronic forward amplitudes fhh are dominated 
by their imaginary parts (see |34j|), we can use the optical theorem to replace 

if hh « -2Mua hN . (23) 

We finally obtain the following expression for the hadron-nucleus cross section: 



.4! 



&hA 



n-1 



[CThNj 



xRe 



//'OC /'OC POO 

d 2 b / dzi / dz 2 ■ ■ ■ / dz„ 
J —oo J Zl J Z„—l 



xp n (b;z 1 • • .^„)exp i 



. Zi — z n 



A 



fi + ^-ir-^K^- 1 ) 

V n=2 I 



(24) 



where 



2 n - 1 (A-n)\ 



Re 



xp n (b; zi...2„)exp i 



d 2 6 / dzi / d^ 2 

^1 — z, n 



dz n 



Zn-l 



A 



(25) 



Note that the exponential factor in Eq. (|24| ) oscillates rapidly if the coherence 
length A of the hadronic scatterer is small. In that case all terms in the series with 
n > 1 approximately vanish and one finds a ha ~ AchN- in the small- x region, 
however, A increases and higher order terms contribute, leading to a reduction of 

&hA- 
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Let us take a closer look at the multiple scattering series (p4"|). For n < i 
we may neglect the recoil motion of the A — n noninteracting nucleons. In the 
absence of nuclear correlations the n-particle density is then approximated by: 



p n (b;2:i ...z n ) 



(26) 



where p is the nuclear one-body density, normalized as / d 3 xp(x) = A. For light 
nuclei only single and double scattering contributions are of importance. The 
above approximation may already be applied for A > 6. Furthermore the validity 
of (^6|) improves with increasing A, since the number of rescatterings n e g- that 
add significantly to OhA grows at most as the nuclear diameter, i.e. n e ff ~ A 1 / 3 . 
For illustration, consider the multiple scattering series fl24|) with a Gaussian 

density p of radius (r 2 ) 1/2 = J%/2 aA 1/3 : 



P(r) 



1 



7r 3/2 fl 3 



cxp 



a 2 A 2 l 3 



We obtain: 



Aa 



hN 



Al/3 a hN A - I 

A * - — - — ; — exp 



2^2/3^ 



.4 



2A 2 



+ 



(27) 



(2? 



We observe that the double scattering contribution adds to the single scattering 
term a negative correction, the magnitude of which grows as A 1 / 3 . Furthermore 
we notice that the shadowing correction decreases rapidly if the coherence length 
of the scatterer becomes small, A < (r 2 ) 1//2 . 



4.2 Shadowing in intermediate-mass and heavy nuclei 

In the previous section we have prepared the tools to calculate total hadron- 
nucleus cross sections <JhA from the respective hadron-nucleon cross sections ct^n- 
We can now proceed to calculate the nuclear structure function as given by 
Eq. (0). 

We will first discuss heavier nuclei, making use of the approximation in 
Eq. (HH), i.e. replacing the n-particle density p n by a product of one-body den- 
sities. We use two 'extreme' parametrizations for these nuclear matter densities: 
a Gaussian shape as in Eq. (|27|) for small A and a square well shape for heavier 
nuclei. In both cases we fit the mean square radii of these density distributions 
to empirical nuclear radii [35] . Note that in earlier calculations [^(J we have used 
realistic densities and included two-nucleon correlations, but we found the result- 
ing corrections in both cases to be systematically very small. With this as an 
input, we can now calculate the shadowing ratios 



R(x,Q 2 



F 2 A ( 



x,Q 2 



AF 2 N (x,Q 2 Y 



(29) 
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Figures ^ and |5] show our calculated ratios for various nuclei, together with exper- 
imental results obtained by the NMC at CERN [|],|| and the E-665 collaboration 
at FNAL who have performed muon scattering measurements focusing on 
the small-x region. We see that for x < 0.1 the ratio fl29|) is generally below 
one, i.e. shadowing occurs. In this x range we can apply the physical picture 
introduced in Section |2|: The virtual photon interacts with the target nucleus 
through hadronic fluctuations. For small x the coherence length A of the fluc- 
tuations becomes large enough to make multiple scattering processes contribute 
significantly. 

However the shadowing ratio (|29|) is not just a function of x but also depends 
(weakly) on Q 2 . We recall from our previous discussion that the value of Q 2 
basically selects that part of the hadron mass spectrum which dominates the in- 
teraction, and hence determines which cross sections ct^jv contribute significantly 
to the multiple scattering series. From Section |3] we note that cihN depends not 
only on the mass // of the qq pair, but also on the distribution of momenta within 
that pair. While the averaged interaction cross sections decrease as log(/z 2 )/^i 2 
with increasing mass, pairs which are asymmetric in the qq phase space interact 
with large cross sections, even for large /i, and therefore produce strong shadow- 
ing. This is the reason for the very weak overall Q 2 dependence of the shadowing 
effect. 

The relevant experiments all operate on fixed targets within a limited range 
of muon energies, hence Q 2 is not an independent parameter but depends on the 
x-range considered. We have taken this dependence into account by inserting 
into our calculation the mean Q 2 values reported for the different x-bins of the 
experiments. With decreasing x the accessible values for Q 2 also become small 
(e.g. at x = 0.005, Q 2 ~ 1 GeV 2 for the NMC experiment from ref. Therefore 
the contributions of the low mass vector mesons p, uj and (f) dominate the observed 
shadowing at x < 0.01 as indicated in Figure 4 and 5. 

The NMC |L],|J has analyzed the Q 2 dependence of shadowing by performing 
linear fits R(x, Q 2 ) ~ a + b In Q 2 to the data for every x-bin. Fig. |G] shows the 
slopes b so obtained in comparison with our calculations. We see that the NMC 
data are compatible with basically no Q 2 dependence. Our calculations give a 
very small positive slope, i.e. a slow decrease in shadowing with increasing Q 2 
which is within the range of the NMC data. 

An NMC analysis of the structure function ratio Sn/C is underway. It com- 
bines data taken at several different muon energies and provides considerably 
better statistics. Figure [7| shows our predictions for this ratio in about the kine- 
matic region to be covered. 

Both the E-665 data on Xenon || and the recent NMC data for Carbon and 
Lithium Q extend to rather small values of x (x < 10~ 3 ). In this region a 
saturation of the shadowing effect becomes apparent, with the ratio eventually 
approaching the 'photon point' i.e. the value observed in the scattering of real 
photons on nuclei (see Section [5] below). Figures |8| a) and b) display the shadowing 
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ratio for Xenon, computed at various fixed values of the energy transfer u, as a 
function of x and a function of Q 2 , respectively. While the onset of shadowing 
is controlled by the coherence length A, which enters as a function of x, one sees 
that the relevant variable for the saturation is Q 2 . As we have already argued, 
variation of Q 2 basically scans the hadronic mass spectrum of the photon. Due to 
the experimental constraints mentioned above, small x in practice always implies 
small Q 2 . Saturation occurs at values of Q 2 less than m 2 , where the interaction 
is dominated by multiple scattering of the p meson. 

Here a remark is in order about contributions from inelastic intermediate 
states to the multiple scattering series (which we have neglected). These are 
significant only at very small x < 10~ 3 and turn out to be small for heavy nuclei. 
Their major contribution to shadowing increases logarithmically with the energy 
transfer v but is independent of Q 2 at small Q 2 < 1 GeV 2 (see Section 4.3 and 
ref. |3B|]). Although the saturation value of R(x,Q 2 ) at x <C 0.1 may therefore 



depend on z/, the onset of the saturation is still controlled by Q 2 . 

Figure |] shows the shadowing ratio YIF^ j AF^ for different nuclei plotted 
against log A at several values of x, together with preliminary NMC data ||37|| . 
One sees that the dependence on A is much weaker than the behavior ~ A 1 / 3 one 
would derive by only considering the double scattering term in (|28|). In fact for 
heavier nuclei, higher order contributions in the multiple scattering series become 
important and partly cancel the effect of the double scattering term. The result 
is a much less pronounced A dependence that can be fitted by the expression 

12F A /AF? « a x + b x In A . (30) 

Figure |10| displays the slopes b x in the different x bins resulting from our calcu- 
lation and those extracted in the preliminary NMC analysis. 

Note, however, that the behavior according to Eq. ( j30"[) cannot be correct in 
the limit of large A. One should rather expect saturation of the shadowing in 
this region, a hint of which can be seen in our results. 

On the whole, our model is able to reproduce the shadowing phenomena in 
deep-inelastic scattering on heavy nuclei remarkably well. With respect to the 
Q 2 dependence, the pending release of the NMC data on Sn may be interesting. 

4.3 Shadowing effects in deuterium 

Shadowing also occurs in deuterium, the most weakly bound nucleus. Although 
small, this effect is of special interest since deep-inelastic scattering from deu- 
terium is used to determine the structure function of the neutron. With the 
assumption of isospin symmetry, the proton and neutron structure functions to- 
gether reveal information on the quark flavor structure of the nucleon. These 
reasons and recent high precision measurements of proton and deuteron struc- 
ture functions and their ratio 0,|3],[8] inspired a lively activity on this topic (see 
e.g. [0,|§). 
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Following our previous discussions, we now calculate the deuteron structure 
function F 2 at small values of x, taking shadowing corrections explicitly into 
account. Consequences for the experimental extraction of the neutron structure 
function F% will then be outlined briefly. 

To calculate F 2 D at x < 0.1 we again need to know the interaction cross 
section <jhD for the scattering of a hadronic fluctuation from the deuteron target 
(see Eq. ([15|)). In addition to incoherent scattering from the two nucleons, o^d 
includes a coherent double scattering correction: 

OTiD = 2a hN - 5a hD . (31) 

From the multiple scattering series in Eqs. (^,|l^) we find 

^^S*"' 1 '' (32) 

In contrast to our treatment of multiple scattering corrections in heavy nuclei in 
Section |4~2| we now take inelastic intermediate states explicitly into account. The 
transition amplitude fhx describes the interaction of the hadronic state h with 
a nucleon by which h is converted into a state X with mass mx- The coherence 
length \(mx) of the hadronic state X is defined as in Eq. (H); we explicitly 
note its dependence on the mass of the propagating state. This coherence length 
enters via the longitudinal form factor Fl of the deuteron, which can be written 
in terms of the deuteron wave function as follows: 

Fl(1/X) 

with ^(r) = ipj=i >M j{r) 

Here u(r) and w(r) are the S- and .D-wave components of the deuteron wave 
function and Xi,Mj * s ^ s s P m P ar t- S 12(f) = 3(oi • f) (cr 2 • f) — (&1 • 02) with 
f = r/|r| is the tensor operator. After taking the average over the target spin we 
find 




1 r°° dz 



2n Jo z 



F L (1/X) = — I — u\z) + w\z) cos ( — )■ (35) 



z 



A 



Our expression for the double scattering correction ( |32|) can be split into an elastic 
(X = h) and an inelastic contribution (X ^ h). As in Section |4J], we assume 
that the amplitudes fhx are strongly peaked in forward direction and dominated 
by their imaginary parts. We may then identify the inelastic contribution with 
the cross section for inelastic diffractive dissociation in the forward direction, 
h + N -> X + N, and obtain: 



1 f d^o"^ nc ^' 

Sa hD = -a 2 hN F L (l/\(m h = fi)) + 8tt J dm\ ' — 



2 u u h^X 



F L [1/X(m x )), (36) 

t=o 
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where t is the squared momentum transfer. 

A well known feature of diffractive dissociation of hadrons and photons is the 
l/m 2 x mass spectrum at large mx (see e.g. ||39||). For hadron h this reads 

( for m\ > n 2 , (37) 



1 d a ofe 



(jfiN dm 2 x dt 



x 

2 



=o m x 



where the constant C ~ 0.1 GeV can be extracted from the experimental analy- 
sis in refs. |40| . We will use Eq. (^) to estimate the inelastic contributions to the 



double scattering correction in Eq. (|36"D , assuming that the l/m x behavior of the 
diffractive cross section sets in at m 2 x = (/i + A) 2 . In high energy hadron-nucleon 
scattering experiments |4(| one finds A typically to be of the order of A ~ 1 GeV. 

With Eqs. flTo] , |3~l"| ) we may now calculate the deuteron structure function F 2 D . 
We use a sample of different deuteron wave functions for this purpose: those 
obtained from the realistic Paris [fyj and Bonn [42| nucleon-nucleon potentials, 
but also - just for comparison - the simple but unrealistic Hulthen ansatz [43]. 
We discuss our results for F 2 D as above in terms of the structure function ratio 

tdDi^ n 2\ _ F?S?h. — 1 _ dFj'jx, Q 2 ) . . 

K [X > Q } ~ 2F 2 N (x, Q 2 ) ~ 2F 2 N (x, Q 2 ) ' m 

In Fig. [TT] we display R D (x, Q 2 ) as a function of x for different values of the 
momentum transfer Q 2 . We observe that R D (x,Q 2 ) < 1 in the range x < 0.1, 
i.e. the characteristic shadowing behavior. The magnitude of the effect is small 
but depends on the deuteron wave function used as an input. For example at 
x = 0.01 and Q 2 = 4 GeV 2 the calculated shadowing effect varies between (1-2)% 
(it amounts to 4% for the naive Hulthen function). 

This sensitivity is a consequence of significant differences, for different po- 
tentials, in the short distance behavior of the deuteron density p(r) = (u 2 (r) + 
w 2 (r))/(47rr 2 ), which determines the longitudinal form factor Fl in Eq. (|35|). In 
Fig. [12] we present p(r) for the various deuteron wave functions, with the densities 
differing considerably for r < 1 fm. On the other hand we note that the region 
r < lfm strongly influences Fl for A > 2fm as can be seen from Eq. fl3~5|) . As 
we have learned in |4.1| , such values of the propagation length control the nuclear 
shadowing effect. 

In Fig. |13| we compare the calculated shadowing effect with and without contri- 
butions from inelastic intermediate states at fixed Q 2 = 4 GeV 2 for the Paris wave 
function. We observe that inelastic states are important only for x < 5 x 10 -3 . For 
example at x = 10 -3 they account for about ~ 20% of the total shadowing effect. 
Their contribution decreases logarithmically with increasing x or, equivalently, 
with decreasing photon energy v. 

We may now briefly justify our statement in section 11 that contributions 
from inelastic intermediate states in multiple scattering are small in heavier nu- 
clei. In analogy to Eq. (|32"D these contributions are proportional to the longi- 
tudinal nuclear form factor of the nucleus which for the example of a Gaussian 
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density reads [BC 



Fi ( 1 / A M~-p(- 5 A&)- (39) 

From Eq. (^) we know that the coherence lengths of intermediate states decrease 
with the invariant mass of the propagating states. Since the invariant mass of the 
hadronic projectile h is always smaller than the mass of the diffractively excited 
inelastic intermediate states, the elastic contribution will naturally dominate the 
multiple scattering process. This dominance is more pronounced as the radius 
(r 2 ) 1,/2 of the nuclear target increases. 

What are the consequences of the shadowing effect in deuterium? As already 
mentioned, the neutron structure function F% is usually extracted from a com- 
parison of the deuteron and the proton structure function. Such an analysis has 
been performed recently by the NMC experiment J3] which investigated the kine- 
matic region x < 0.1 with high accuracy. In this analysis, however, effects from 
nuclear shadowing in deuterium have been ignored; the difference of proton and 
neutron structure function was obtained by simply taking 

-^2 P nmc = ( F 2 - ^")nmc = 2-Ff - F 2 D . (40) 

However, for small x shadowing must be taken into account. The true structure 
function should therefore read: 

F V-n = F P _ F n = 2F p _ {p D + 5F D } = F| -n Mc _ (41) 

The shadowing correction 5F® reduces the result with respect to the quoted dif- 
ference nmc ( an d correspondingly the true neutron structure function should 
be larger than the value obtained by the NMC). 

The full symbols in Fig. IT4] display the original NMC data for the difference 



nmc as we U as corrected results with the shadowing term subtracted for small 
x. In our calculation we used Q 2 = 4GeV 2 , as the NMC analysis operates with 
structure functions interpolated to this value, and again a set of different deuteron 
wave functions. We notice that for small x the structure function difference 
becomes small, so that the relative size of the correction is of the order of 100%. 
This is in good agreement with the expectation that any deviation of i^NMC 
from zero in this x region should be mostly due to the shadowing effect. 

Let us next consider the integral over the difference of proton and neutron 
structure functions: 

I G (x, 1) = / — (FJV) - F 2 V)) • (42) 

J x X 

In the parton model Iq can be rewritten in terms of quark distributions: 

I G (x, 1) = ~ f 1 dx' (u v (x') - d v (x')) + lJ l dx' (u(x') - d(x')) . (43) 

6 Jx O Jx v ' 
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The up and down valence quark distributions u v = u — u and d v = d — d in 
the proton are given by the difference of the respective quark and antiquark 
distributions. To obtain Eq. (f43|) we have used isospin symmetry to relate proton 
and neutron quark distributions. 

If the first moments of the up and down sea quark distributions are approxi- 
mately equal - this is fulfilled trivially if one assumes the sea to be SU(2)-flavor 
symmetric - the second term in Eq. (|43|) vanishes for x — > 0, and one arrives at 



the Gottfried sum rule [44 



I G = I G (0, l) = \f l dx' K(x') - d v (x')) = 1 (44) 

6 JO 3 

Without taking shadowing into account, the NMC found 7^ xp '(0,l) = 0.235 ± 
0.026 ||. This includes contributions from the unmeasured regions x > 0.8 
and x < 0.004. A smooth extrapolation of -F^/Fj' for x — > 1 yields Jg(0.8, 1) = 
0.001±0.001, while within conventional Regge theory J G (0, 0.004) = 0.013±0.005 
is found J|. 

The deviation of the Gottfried sum I G xp ' from the naive expectation 1/3 has 



been the target of some activity, see e.g. f45[. 

Let us now consider the impact of our corrections for deuteron shadowing 
on the extraction of the Gottfried sum. As discussed above, for x < 0.1 they 
significantly reduce the difference i 7 f~ n with respect to the values used by the 
NMC collaboration. 



The open symbols in Fig. |TJ] show the NMC results for the integrals 
Ig{%, 0.8) together with our corrected values. In Tab. ^| we give the corrections 
A7g(0.004, 0.1) due to shadowing obtained for different deuteron wave functions. 
They reduce the extracted value for the Gottfried sum by 10% or more. This 
sizable correction is due to the fact that the structure function difference in the 
integral in (|42| ) is weighted by a factor 1/x. We see that the shadowing correc- 
tions further enhance the deviation of the Gottfried sum from the naive value 
I G = 1/3. 

The E665 collaboration has also measured the structure function ratio F^/F^, 
for which preliminary data are now available ||. Their x-range extends down to 
10~ 6 , with the average Q 2 strongly dependent on x ((Q 2 ) = 0.002 GeV 2 for the 
lowest x bin 10 -6 < x < 10~ 5 ). Figure [15| shows the ratio F2/2F2 as obtained 
by the E665 group. The small-x behavior of these data is in excellent agreement 



with the results of our model calculation. In comparing with Figs. |TT] and [13 
we note that the logarithmic growth of shadowing for x — ■> displayed there 
is a result of plotting the ratio at constant Q 2 , which implies a photon energy 
~ 1/x. In this experiment the muon (and hence the photon) energy is limited 
and therefore the contributions from inelastic intermediate states are bounded. 
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5 Nuclear shadowing of real photons at high energy 



In deep-inelastic scattering, the lepton beam is a source of highly energetic virtual 
photons. We have conveniently written the structure function in terms of a 
total cross section for the interaction of the virtual photon with the target (see 
Eq. D. In the limit Q 2 -> this cross section becomes physical, describing real 
photoproduction processes. 

While hadronic vacuum fluctuations decouple from the photon in the limit 
Q 2 — > 0, so as to keep the real photon massless, shadowing still occurs via the 
production of quark-antiquark pairs on nucleons and their subsequent propaga- 
tion through the nucleus. The coherence length of such hadronic states of mass 
/i, produced by a real photon of energy u, is 

A = % (46) 

fl 2 

This becomes largest for those states with the smallest mass. For the p meson 
with = m 2 we find that A reaches typical internucleon distances for v > 3 GeV. 
We therefore expect multiple scattering to reduce the cross section for nuclear 
photoproduction o 1 a with respect to A times the photoproduction cross section 
on free nucleons a lN . 

Let us discuss this more quantitatively. The expression for a 1 N at large v is 
obtained from Eq. ([13]) taking the limit Q 2 — > 0. This yields 



47ra em f°° f 1 n cont, (u 2 ) 
V — ^-a VN + 47ro em / dfi 2 / da a hN (fi 2 ,a). (46) 

v 9v M Jo V 



The VMD term (~ 90/ib) is about twice as large as the continuum contribution. 

As in the case of deep-inelastic scattering, we can calculate the cross section 
for photoproduction on nuclei by replacing the hadron-nucleon cross sections in 
Eq. (f|6|) by the corresponding hadron-nucleus cross sections, obtained via the 
multiple scattering formalism as outlined in Q| . To discuss nuclear effects it is 
common to consider the ratio 

A Aa lN 

In our calculation of a 1 N we use parameters which were fixed through our fit to the 
nucleon structure function in Section |3|. To calculate <t 7j 4 we employ similar 
nuclear densities as in the case of nuclear deep-inelastic scattering in Sec. |4.2| . 

The A-dependence of photoproduction on nuclei has been measured by sev- 
eral groups [^]. Figure [16] displays our results of the shadowing ratio A eS /A 
calculated for C, Cu and Pb target nuclei, together with various experimental 
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data. We observe significant shadowing for u > 2 GeV. The effect is well de- 
scribed within our model. Its results are quite similar to those of earlier VMD 



calculations [13], which should not be much of a surprise, due to the dominant 



role of the VMD term noted above. 

Shadowing grows stronger for higher photon energies v, eventually approach- 
ing some saturation value, apart from logarithmic corrections due to inelastic 
contributions to multiple scattering (see discussion in Sections 4.2 and 4.3). This 
is due to the fact that the coherence length A governing multiple scattering pro- 
cesses is now directly proportional to v. 

The validity of the picture developed here is restricted to v > 2 GeV. For 
smaller photon energies, photonuclear dynamics is governed by the excitation 
and propagation of nucleon resonances in nuclei. 



6 Summary 

Shadowing at small values x < 0.1 of the Bjorken variable is the most prominent 
nuclear effect seen in deep-inelastic lepton scattering from nuclear targets. We 
have developed a phenomenology of nuclear shadowing, expressed in the labo- 
ratory frame, which makes use of the full hadronic spectrum of virtual photons. 
Our framework unifies the vector meson dominance picture with the concept of 
color transparency applied to the quark-antiquark continuum part of the photon 
spectral function. 

Our results are summarized as follows: 

i) In our lab frame approach shadowing arises from the coherent multiple 
scattering of quark-antiquark fluctuations through the nuclear target. A 
satisfactory description of nuclear structure functions at small x is achieved 
for a large variety of nuclei, both light and heavy. 

ii) Vector mesons dominate at small Q 2 < 1 GeV 2 . At large Q 2 the quark- 
antiquark continuum becomes important. The combination of both reso- 
nant and continuum parts of the hadronic photon spectrum is crucial in 
order to obtain the almost negligible overall Q 2 dependence of the shadow- 
ing effect, whereas vector meson dominance alone would imply decreasing 
shadowing with increasing Q 2 . 

iii) Contributions to shadowing from inelastic intermediate states in the mul- 
tiple scattering chain give only small corrections for heavy nuclei, but they 
need to be taken into account for light nuclei, in particular for the deuteron. 

iv) Consistency is found with the observed shadowing for interactions of real 
photons with nuclei at high energies. 
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v) Special emphasis has been directed to shadowing effects in deuterium. In 
the extraction of the neutron structure function from deuteron data at 
small x, such effects must be taken into account carefully. We find strong 
sensitivity to the short distance behavior of the deuteron wave function; the 
use of realistic nucleon-nucleon potentials is therefore absolutely necessary 
for a reliable description. Shadowing effects of 1-2% imply corrections of 
the order of 10% in the Gottfried sum. This correction further increases 
the already established discrepancy with the naive parton model. 

We would like to thank S. Kulagin, W. Melnitchouk, N. N. Nikolaev and A. W. Thomas 
for helpful discussions and comments. 
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Tables 



V m v I MeV |§ g v /4n a VN / mb 



p 769.9±0.8 2.0 

u 781.9±0.1 23.1 

<p 1019.41±0.01 13.2 

J/ip 3096.88±0.04 10.5 

ip' 3686.0±0.1 30.6 



22-27 
25-27 
9-12 
2.2 ±0.7 P7fl 



HI 
El 







1.3 



Table 1: Vector meson properties: masses, couplings to the photon, and total 
vector meson-nucleon cross sections. The coupling constants gy are derived from 
the V —>■ e + e~ decay widths 



AI G (0.004, 0.1) 
I G xp -(0.004, 0.8) [1 Bonn(l) Bonn(2) Paris 
0.221 ±0.021 -0.022 -0.039 -0.017 



Table 2: Shadowing corrections AIg for the Gottfried sum obtained for vari- 
ous deuteron wave functions compared to the experimental value '. Here 
"Bonn(l)" refers to the one-boson exchange Bonn potential, "Bonn(2)" is the 
full potential including explicit two-pion exchange etc. |4"2"|] . The Paris potential 
is taken from ref. ||41|| . 
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Figure captions 



1. The two possible time orderings for the interaction of a (virtual) photon 
with a nucleon or nuclear target: (a) the photon hits a quark in the target, 
(b) the photon creates a qq pair that subsequently interacts with the target. 



2. Nucleon structure function for small x plotted against Q 2 . The solid line 
is the full result of our calculation. The contribution of vector mesons is 
indicated by the dashed line. We compare to NMC data from ref. ||. 

3. Contribution to the multiple scattering series: the hadronic projectile 
scatters from n nucleons inside the target nucleus. 

4. Our results for shadowing in He, Li, C, and Ca compared to available 
experimental data [I],f3],[7]. The dashed curves show the shadowing caused 
by the vector mesons p, u and only. 

5. Shadowing in Xenon. Data are from the FNAL E-665 experiment [^, ||. 
The vector meson contribution is shown by the dashed curve. 

6. The slope b = d-R/dlnQ 2 indicating the Q 2 dependence of the shadowing 
ratio for He, Li, C, and Ca extracted by the NMC fl],§] for various x-bins 
together with our results. 

7. Q 2 dependence of the shadowing ratio Sn/C as predicted by our model, in 
the region to be covered by recent NMC data. 

8. The shadowing ratio in Xe (a) as a function of x at fixed v\ (b) as a function 
of Q 2 at fixed v. 

9. The shadowing ratio as a function of the nuclear mass number A for several 



x-bins. Experimental data are preliminary NMC results from [37 



10. Our results for the slopes b x of the A dependence relative to Carbon ob- 
tained from Eq. ( |30| ) compared to the NMC data in ref. |37 . 



11. Shadowing ratios in deuterium plotted against x at different values of Q 2 
and for different deuteron wave functions. Here "Bonn(l)" refers to the one- 
boson exchange Bonn potential, "Bonn(2)" is the full potential including 
explicit two-pion exchange etc. f42|| . The Paris potential is taken from 
ref. gTJ. 



12. The density p(r) = (u 2 (r) +w 2 (r))/4nr 2 corresponding to different param- 
etrizations of the deuteron wave function. 
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13. The shadowing ratio for deuterium at fixed Q 2 = 4 GeV 2 calculated for the 
Paris wave function. The full line includes inelastic intermediate states. 
The dashed curve was obtained by taking only elastic intermediate states 
into account. 

14. Gottfried sum and shadowing: the filled symbols display the structure func- 
tion difference -F 2 p_n (:r), with the circles representing the original NMC 
data ||, squares and diamonds include shadowing correction using differ- 
ent deuteron wave functions. Open symbols show the respective values for 
the integral Iq(x, 0.8). 

15. Our result for R D (x,Q 2 ) compared to recent FNAL E-665 || and NMC 
data @] for the ratio F$j2F%. 

16. The shadowing ratio for the absorption of real photons on nuclei as calcu- 



lated in our model and measured by various groups |46| plotted against the 
photon energy v. 
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